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Context. Cold fronts have been observed in a large number of galaxy clusters. Understanding their nature and origin is of primary 
importance for the investigation of the internal dynamics of clusters. 

Aims. To gain insight on the nature of these features, we carry out a statistical investigation of their occurrence in a sample of galaxy 
clusters observed with XMM-Newton and we correlate their presence with different cluster properties. 

Methods. We have selected a sample of 45 clusters starting from the B55 flux limited sample by Edge et al. (1990) and performed a 
systematic search of cold fronts. 

Results. We find that a large fraction of clusters host at least one cold front. Cold fronts are easily detected in all systems that are 
manifestly undergoing a merger event in the plane of the sky while the presence of such features in the remaining clusters is related to 
the presence of a steep entropy gradient, in agreement with theoretical expectations. Assuming that cold fronts in cool core clusters 
are triggered by minor merger events, we estimate a minimum of 1/3 merging events per halo per Gyr. 
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1. Introduction 

The unprecedented angular resolution of the X-ray telescope 
Chandra led to the discovery of several new phenomena within 
various astrophysical systems. One of these are cold fronts de- 
tected in galaxy clusters. Initially ob served in merging clus - 
ters, the prototypes are found in A2142 ( Markevitch et al. 2000), 
A3667 ( Vikhli njneFal.11200 1 alibi: IVikhUnin & Markevitch| [2002) 
and 1E0657-56 dMarkevitch et al.ll2002l) . All these systems fea- 
ture very sharp discontinuities in their X-ray images where the 
drop of the surface brightness (and correspondingly of the gas 
density) is accompanied by a jump in the gas temperature, with 
the denser region colder than the more rarefied region, unlike 
shock fronts. For this reason, these features have been dubbed 
"cold fronts" (Vik hlinin et al.ll2001al) . The density and the tem- 
perature discontinuities have similar amplitude so that pressure 
is approximatively continuous across the front. 

Cold fronts have been initially interpreted as the edge of 
the cool core of a merging substructure which has survived 
the merger and is rap idly moving through the ambient gas 
dMarkevitch et all2000l) . 

Cold fronts have successively bee n detected in the cores 
of some relaxed cl usters (e.g. A1795: i Markevitch et al.l l20Qlt 
RX J 1720. 1+2638: Maz zotta et alJ200U A4 96: lDupke & White! 
120031: 2A 0335+096- lMazzoita et alj|2003l) and to date a large 
number of relaxed systems are known to host one. Since the 
presence of cold fronts in cool cores provides evidence of 
gas motions and possibly of departures from hydrostatic equi- 
librium, understanding the nature of such a widespread phe- 
nomenon is mandatory to characterize the dynamics of galaxy 
clusters. High resolution hydrodynamical simulations are, at 
present, the main technique to investigate the mechanisms gen- 
erating cold fronts. Indeed, cold front features could already 
be detect ed in simulations published prior to the launch of 
Chandra dRoettiger et alj I1997L [[998). After cold fronts dis- 



covery, several hydrodynamical simulations have been devel- 
oped to model the effect of the ram-pressure stripping in a 
merger event and the formation of the cold front feature in 
merging clusters dHeinz et al.l 120031: iNagai & Kravtsovl [2003; 
Math is et al.l [2005). Several simulations have also been em- 
ployed to understand the origin of cold fronts in relaxed non- 
merg i ng clusters (e.g.lChurazoy et al. 2003; Tittle v & Henriksenl 



Tglt 

2005; As casibar & Markevitchll2006l) . 
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The emerging picture (jAscasibar & Marke vitchl 120061 : see 
also IMarkevitch & Vikhlinin 2007 for a review) is that cold 
fronts arise during major merging events through ram-pressure 
stripping mechanisms which induce the discontinuity among the 
merging dense subcluster and the less dense surrounding ICM. 
In relaxed clusters, the cold fronts features are induced by mi- 
nor merger events which produce a disturbance on the gas in the 
core, displace it from the center of the potential well and de- 
couple it from the underlying dark matter through ram-pressure. 
Subsequently, a sloshing mechanism sets in, generating cold 
fronts. The necessary condition for triggering this mechanism is 
the presence of a steep entropy profile for the central gas which 
is generally fulfilled at the center of relaxed cool core clusters. 

Cold fronts are at present observed in a large number of 
galaxy clusters. [Markevitch e t al. (2003) analyzed a sample of 
37 relaxed clusters observed with Chandra showing that cold 
fronts are pres ent in the ma j ority o f the cores of relaxed clus- 
ters. Recently, lOwers et al.1 d2009t) characterized a sample of 
nine cold fronts with quantitative measurements of the thermo- 
dynamic discontinuities across the edges and associated the pres- 
ence of a cold front with evidence of merger activity. 
While many objects have been studied in detail to understand the 
nature of cold fronts, we still lack a systematic investigation of 
the characteristics of these phenomena and of their host clusters 
through a large sample. The aim of this paper is to perform a 
systematic search of cold fronts in a representative sample and 
to investigate the properties of their parent clusters. Such a study 
is necessary to inspect the nature and origin of cold fronts and 
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eventually to test the reliability of the picture emerging from the 
simulations. 

The sa mple is selected s tarting from the B55 flux limited 
sample by Ed ge et all (Q.990). We use for our analysis XMM- 
Newton data. In spite of its limited spatial resolution with respect 
to Chandra, XMM-Newton has the positive attribute of having a 
large field of view, allowing a significant coverage of most of 
the clusters. In most cases, the clusters are inside the EPIC field 
of view up to a radius J 0.3ri8o, allowing the characterization 
of the main thermodynamical properties well beyond the core 
regions. Additionally, the XMM-Newton large collecting area al- 
lows a good statistics for a large number of objects. 

Among the several physical properties characterizing the in- 
tracluster medium, we focus our attention on entropy. Entropy 
plays a key role in describing the thermodynamical state of 
the ICM, its distribution is a signature of the thermodynam- 
ical history of the cluster and it is also intimately related 
to the non-grav i tation al processes which may have occurred 
(fVoit et al.1120021 l2005h . Moreover, as previously stressed, sim- 
ulations highlight how the steep gas entropy profile is a nec- 
essary condition for the onset of the sloshing mechanism and 
therefore for the presence of c old fronts in cool core clusters 
dAscasibar & Markevitchll2006l) . 

The structure of the paper is the following. In §|2]we describe 
the sample of clusters that we have analyzed and in §[3] we pro- 
vide details about the data reduction. Then we describe (see § [4j 
the algorithm used for the systematic search of cold fronts in the 
cluster sample. We present our results about the occurrence and 
the origin of cold fronts in § [5] and we discuss them in § [6j We 
summarize our findings in §|7] 

We adopt a ACDM cosmology with Q m = 0.3, Q A = 0.7, 
and Ho = 70kms- 1 Mpc 1 . 



2. The sample 

We use as a refere nce starting sample the flux limited sample by 
lEdge et al.ldl990h . It includes 55 objects with flux f x > 1.7-10 -11 
erg cirT 2 s _1 in the 2-10 keV energy band and is 90% complete. 
All the clusters are located at redshift z < 0.2. 

Starting from this sample, we select all the clusters avail- 
able in the XMM-Newton public archive. At the time of writing, 
Ophiuchus was not publicly available and the objects A2244, 
A644 had not been observed with XMM-Newton. For A754 
we use the long observations we obtained from A07 (P.I.: 
A. Leccardi). Observations for 3C129, A2142, A2147, A1736, 
A3391 are highly affected by soft protons: since their final good 
exposure time after cleaning procedures is generally below 5 ks 
for MOS and < 0.5 ks for pn, these clusters have been excised 
from the sample. For clusters having more than one observation, 
we eliminate those observations with high soft protons contam- 
ination. We also exclude Virgo and Coma: their extension does 
not allow a significant coverage within the EPIC field of view. 

Our final sample is reduced to 45 objects. In TableQ]we filed 
the list of clusters belonging to our sample. 

The excision of a number of clusters invalidates the com- 
pleteness of the sample. To verify if the final sample is represen- 
tative of the cluster population, we inspect the distribution of the 
main cluster observables. We built the histograms for redshift, 
X-ray luminosity and temperature (see Fig.Q} both for the start- 
ing sample (light grey) and for the exclud ed clusters (dark grey). 
For X-ray luminosities (Lx) we refer to Reiprich & Bohringer 
(2002), while temperatures are taken from iPeres et alj dl998). 
The histograms show that excluded clusters do not introduce any 
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Fig. 1. Distributions of redshifts ( a), X-ray luminos i ties (b) and 
tempe ratures (c) for the starting (Ed ge et alj|1990t ; IPeres et all 
1998) clusters sample (light grey) and for the excised objects 
(dark grey). The rejection of these objects does not introduce 
any obvious bias on the remaining subsample. 



obvious bias. We conclude that, even if the final adopted sample 
is not complete, it is representative of the cluster population. 



3. Data reduction 

Observation Data Files (ODF) were retrieved from the XMM- 
Newton archive and processed in a standard way with the 
Science Analysis System (SAS) v6.1. 

We apply the standard selection #XMMEA_EM to the MOS 
event list (#XMMEA_EP for pn) to automatically filter out artefact 
events. The soft prot ons cleaning was performed using a double 
filtering process (see Leccardi & Molendi 2008; Pratt & Arnaud 
2002). The adoption of a threshold level and the exclusion of 
light curve intervals above the selected threshold allows the re- 
jection of most flare events. In practice, we extract the light curve 
in the 10 - 12 keV (10-13 keV) energy band for MOS (pn) using 
100 second bins. We apply a threshold of 0.20 cts s _1 for MOS 
and 0.50 cts s for pn to generate the filtered event file. However 
softer flares may exist such that their contribution above 10 keV 
is negligible. To remove t his flare contamin ation, we apply the 
3cr clipping method (see iMartv et alj|2003l) : we extract a his- 
togram of the light curve in the 2-5 keV band and fit this his- 
togram with a Gaussian distribution. Since most flares have al- 
ready been rejected in the previous step, the fit is usually very 
good. We then apply a threshold at the +3cr level and generate 
the filtered event file. After soft proton cleaning, we filter the 
event file according to FLAG (FLAG==0) and PATTERN criteria 
(PATTERN<12). 

To systematically search surface brightness discontinuities, 
we build for each cluster, the EPIC flux map: MOS1 + MOS2 + 
pn (pn images are corrected for out of time events). This flux 
image is computed in the 0.4 - 2 keV band following a pro- 
cedure similar to the one described in iBaldi et al.1 ([2002, see 
also iRossetti et al.1 120071) . We sum up the MOS1, MOS2 and 
pn source images to obtain the total source map, S epic, and we 
compute EPIC source exposure map, expEPic, by summing the 
source exposure maps of each detector. The EPIC count rate im- 
age is then defined as crEPic = S epic I ex PEPic- Count rates are 
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cluster name 


redshift 


cluster name 


redshift 


Centaurus 


0.0114 


A85 


0.0551 


A 1060 


0.0126 


A3532 


0.0554 


A262 


0.0163 


A3667 


0.0556 


AWM7 


0.0172 


A2319 


0.0557 


Perseus 


0.0176 


Cygnus A 


0.0561 


A 1367 


0.0220 


A2256 


0.0581 


A4038 


0.0300 


A3266 


0.0589 


A2199 


0.0302 


A3158 


0.0597 


A496 


0.0329 


A1795 


0.0625 


2A0335 


0.0349 


A399 


0.0718 


A2063 


0.0349 


A2065 


0.0726 


A2052 


0.0355 


A401 


0.0737 


A576 


0.0389 


A3112 


0.0750 


A3571 


0.0391 


A2029 


0.0773 


A119 


0.0442 


A2255 


0.0806 


MKW3s 


0.0450 


A 1650 


0.0838 


A 1644 


0.0473 


A1651 


0.0849 


A4059 


0.0475 


A2597 


0.0852 


A3558 


0.0480 


A478 


0.0881 


A3562 


0.0490 


PKS0745 


0.1028 


Triang. Aus. 


0.0510 


A2204 


0.1523 


Hydra A 


0.0538 


A 1689 


0.1832 


A754 


0.0542 







then converted to flux through the total conversion factor c/epic 
derived following the formula: 

expEPic _ exp MO s \ exp M0S 2 exp p „ 

c fEPIC c/moSI c/moS2 Cf pn 

where cf M os\,cf MOS2 ,cf pn and exp M0 si,exp M 0S2,expp n are 
the conversion factors and the exposure maps of the three instru- 
ments. The EPIC source flux image is obtained using the relation 

fx epic = c/epic • crEPic- 

To remove the quiescent particle induced background and 
the cosmic background component we need EPIC background 
flux images. We use a large collection of background data, such 
as long observations of blank sky fields. We use 9 blank fields 
for a total exposure time of ~ 300 ks) selected by our own group 
Leccardi & Molendill2008l) . We compute the EPIC background 
flux image from the MOS1, MOS2 and pn background images 
and exposure maps using the same method applied to the EPIC 
source flux image. By subtracting the EPIC background flux im- 
age from the source flux image we derive an EPIC net flux im- 
age in units of 10~ 15 erg cm -2 s _1 pixel -1 (one pixel is 8.5 x 8.5 
arcsec 2 ). The net maps are used for the construction of the sur- 
face brightness profiles. 

To search for cold fronts, we also need the temperature maps 
of all the clusters of our sample. We adopt a modified version 
of the "ad aptive binning + B road Band Fitting" algorithm de- 
scribed in Rossetti et al. (2007), where we have substituted the 
Cappellari & Copini l2003l) adaptive binnin g algorithm with its 
modified version bv lDiehl & Statlerl (120061) . 



4. Systematic search of cold fronts 

4. 1 . The detection algorithm 

In this section, we identify a suitable method to detect cold 
fronts, well aware that any general selection criterion will have 
some limitations and can introduce spurious effects, so that some 
cold fronts may be missed and some features may be classified 



as cold fronts although they are not. We will address this point 
in detail in § Roland in §1511 

A cold front is characterized by the presence of a sharp 
decrease in the surface brightness (SB) profile typically ac- 
companied by a rise of the gas temperature. We initially per- 
form a systematic search of surface brightness discontinuities in 
the cluster sample and generate a list of candidate cold fronts. 
Subsequently, we examine the gas temperature behavior across 
the discontinuity, in order to rule out the hypothesis that the de- 
tected discontinuity is a shock front. 

We developed an algorithm to perform the systematic search 
of surface brightness discontinuities. We start from the EPIC flux 
maps (see § O that we have built for our sample and we divide 
each cluster map in 30° wide sectors centered on the SB peak. 
In most cases, we detect SB discontinuities in several consec- 
utive sectors, suggesting that most cold fronts have an angular 
extension significantly larger than 30°. Consequently, unless the 
statistics is particularly low, fixed angular ranges can be used to 
find discontinuities. A possible bias on the detectability of cold 
fronts due to the width of the angular sectors will be discussed 
in § El For clusters having a low statistics or where the possible 
cold front is located near the cluster center (e.g. A262, A1795, 
A2199), we use an "ad hoc" choice of the sectors (45° wide or 
even larger) to reveal the discontinuity in the surface brightness 
profiles. We build the surface brightness profile for each sector 
using the cluster X-ray emission peak as center. Sometimes, for 
merging or irregular clusters, a visual inspection of the images 
may suggest a different center, better suited to detect a sharp 
decrease of the surface brightness. In Fig. [2] we show the sur- 
face brightness images of A2319 and A1367 as an example. The 
black circles mark the centers we adopted to build the profiles. 

For each cluster a set of profiles is obtained. While for some 
clusters (e.g. Centaurus, A496, 2A0335+096) the presence of 
a surface brightness discontinuity is apparent, in other clusters 
(e.g. A262) the surface brightness profile is not as sharp (see Fig. 
|3}. Projection effects and resolution limits smooth the profiles: 
the surface brightness discontinuity will appear as a steepening 
of the profile in the radial range around the jump radius. In the 



4 



Simona Ghizzardi et al.: Cold fronts in galaxy clusters 




Fig. 2. Surface brightness map for A2319 (top panel) and A 1367 
(bottom panel). The maps have been smoothed for a better visual 
inspection. The black circles show the position of the centers 
we chose for the radial profiles. They roughly match the centers 
of curvature of the candidate discontinuity (SE-E for A2319 and 
S for A 1367) and do not match the X-ray emission peak. This 
choice allows a better characterization of the jump in the sur- 
face brightness profiles. Full resolution figures are available at: 
http://www.iasf-milano.inaf.it/~simona/pub/coldfronts/ghizzardi 



approximation where profiles are described by power laws, the 
slopes measure the steepness of the profile. We use the power 
law slopes to characterize the surface brightness discontinuities. 
We identify for each cluster the possible discontinuity region 
with a visual inspection of the profile and of the image. We mark 
this region with the letter D and we set SB oc r~ a ° in the cor- 
responding radial range (see upper right panel of Fig. [3] as an 
example). We compare the slope we find in this region with the 
slope obtained fitting the profile with the power law SB oc r~ am 
in the nearby (inner and/or outer) region or in other sectors of the 
cluster where no irregularities in the surface brightness profiles 
are present. The difference of the slopes Aa - ao - ond quan- 
tifies the steepening of the profile. We require that Aa > 0.4 
to classify a region as discontinuous. The choice of this thresh- 
old relies on phenomenological considerations. All the jumps 



we measured have Aa values well above 0.5, while for regions 
without discontinuities Aa values are below 0.2. 

Examples of different surface brightness profiles are reported 
in Fig. [3] Centaurus cluster profiles (top left panel in the fig- 
ure) steepen significantly between 50" - 100" in the NE sector 
and between 170" - 210" in the W-NW sector. In A262 and in 
A 1060 (top right panel and bottom left panel respectively) the 
discontinuities are not as apparent as in Centaurus. Finally, the 
profile for AWM7 (bottom right panel) does not show any irreg- 
ularity. We fit all the profiles in the different radial ranges with 
power laws. In Table [2] we report, for each region of the four 
clusters shown in Fig. [3] the ranges used for the fits, the slopes 
of the best fits, and the associated Aa values. 

For the Centaurus cluster, we consider two interesting sec- 
tors, SE and W-NW of the cluster core. The slopes we find for 
the W-NW sector (specifically, 120° - 150°, where the angles 
are measured in an anticlockwise direction from East) are 0.41, 
2.59 and 0.84 for the radial ranges: 10" - 50", 50" - 100", 100" 

- 240" respectively. The Aa for the central region is 2.19 with 
respect to the innermost region and 1 .75 with respect to the outer 
one. This is obviously classified as a discontinuity and is a can- 
didate cold front. A similar analysis allows to assess that there is 
a discontinuity in the 170" - 210" radial range in the 30° - 60° 
sector. The quality of the Centaurus cluster profiles is extremely 
high thanks to its proximity and to the long observations (170 
ks in the public archive at the time of writing) so that, even at 
the margin of the cold front located west of the core, from -60° 
to 60°, the discontinuity is still visible in the profile. The Aa we 
find for this case, 1.60, is smaller than in the previous case, but 
still high. 

In the top right panel, we show the profile for the SW-W 
region (-45° - 0°) in A262, where a discontinuity at ~ 60", al- 
beit not very sharp, is identifiable (Aa is slightly smaller than 
1). In A1060 (bottom left panel) the putative discontinuity is 
around 2'. As shown in Table [2] the analysis of the slopes pro- 
vide Aa = 0.02. As a consequence, this feature is not classified 
as discontinuity. However, a different result could be obtained 
with a slightly different choice of the radial ranges used for the 
fits. If we restrict the range of the discontinuity region to 110" 

- 130" (3 points for the fit) and we choose 150" - 250" for the 
outer region (disregarding the 4 points immediately after the dis- 
continuity where the profile is flat) Aa increases to 0.55 and this 
feature could be considered as a discontinuity. We believe that 
this choice of the radial range is rather extreme. In addition, we 
note that the temperature profile shows no variations in the same 

p^jgion. As a general rule, if Aa satisfies the discontinuity condi- 
tion only for an ad hoc choice of the radial range, we exclude it 
from the list of candidate cold fronts. Finally, in the bottom-right 
panel in Fig. [3] AWM7 shows a regular behavior. In the figure 
we report one region, but AWM7 is regular in all its sectors and 
all profiles are similar; for this cluster we find no discontinuities. 

The procedure described thus far detects surface brightness dis- 
continuities and provides a list of candidate cold fronts. To up- 
grade a discontinuity to a cold front, we need to verify the be- 
havior of the temperature profile across the surface brightness 
jump. T o this aim, we build the binned temperature maps (see 
§ l3l and iRossetti et alJ 120071) for all the clusters of our sample. 
From these maps we derive the temperature profiles plotting all 
the bins whose barycentre is inside the sector hosting the surface 
brightness discontinuity. In none of the candidate cold fronts we 
observe a sharp decrease in the temperature profile as would be 
expected for a shock front. Almost all the surface brightness dis- 
continuities that we find feature a sharp gas temperature rise. In 
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Fig. 3. Surface brightness profiles for some sectors of four clusters of our sample. For the Centaurus cluster we plot two interesting 
sectors: SE (filled circles) and N-NW (open diamonds) of the cluster core. The figure shows that profiles may have different behav- 
iors. While discontinuities are apparent in some clusters (e.g. Centaurus cluster), in others they are not as sharp. Some systems (e.g. 
AWM7) have a regular profile. In upper right panel (A262) we plot the ranges used to fit the profile with power laws. The flag D 
marks the discontinuity region and ND marks the adjacent (inner and outer) regions. In all the panels the solid lines represent power 
law best fits (see text and Table[2]for details). In these plots, surface brightness is given in 1CT 15 erg cirT 2 s arcsec~ 2 units. 



some cases the temperature rises smoothly with no jumps. We 
remark that the cold front feature does not necessary require a 
temperature jump, since the thermal pressure of the gas inside 
the cold front is balanced by the sum of the thermal and ram 
pressures outside. This can be achieved also with a slow contin- 
uous rise of the temperature across the discontinuity. 

4.2. Notes on individual clusters 

In Tab. [3] we list the clusters hosting one or more cold fronts. 
For each cluster, the table provides the center used to build the 
SB profiles, the cold front, the radial and azimuthal position of 
the cold front and Aa. The Aa reported for each cold front is the 
mean value obtained averaging over the different sectors hosting 
the discontinuities. 

As already remarked, the procedure adopted to classify cold 
fronts can fail in finding the discontinuities or provide some 
doubtful cases. Hence, comments are required for some individ- 
ual systems. 

- A4059: a visual inspection of the surface brightness map of 
A4059 hints to a possible cold front in the SW sector ~ 30" 
from the peak, but the discontinuity is barely visible in the 
profiles and we did not succeed in fitting it with power laws 
and deriving Aa. We consider this as an unclassified case. 

- A85: This object has two possible cold fronts. The former 
lies in the NW sector ~ 80" from the X-ray peak. In this 



sector Aa = 0.75, above our threshold. However, the sector 
of the cold front is very narrow (30°) and no discontinuities 
are detected in the nearby sectors. Table [2] shows that the 
cold fronts widths generally range from 60° to 120°. Such 
a tiny extension for a cold front is unusual and we prefer to 
consider this as an unclassified case. Another cold front is 
present in a small subclump located 8' south of the cluster 
and moving north towards the main structure dDurret et al.l 
1998; Ke mpner et al.l2002l) . This is labeled as A85* in Table 
[3]to distinguish it from the unclassified cold front of the main 
cluster. 

- A2052: the analysis of the discontinuities in A2052 is com- 
plex because of the presence of the bright she lls surrounding 
the X-ray cavities dBlanton et alj|2003l 1200 ll) . Remarkably, 
sharp decrements of the surface brightness profiles are de- 
tected just outside the rims, at distances of about 40"-50" 
(rims are at 10"-30" from the center). It is difficult to es- 
tablish whether such sharp drops are real discontinuities or 
if they are associated to the bright shells. Recent results 
obtain ed from a deep Chandra observation dBlanton et alj 
120071) show that some cavities and ripples are present in this 
cluster, similarl y to what has been observed in Perseus clus - 
ters and M87 dFabian et al.ll2006l, 120031; iForman et aTll2007l) . 
Some weak shocks may also be present. The presence of cold 
fronts is unclear and we consider this cluster as unclassified. 

- Hydra A: in Hydra A, we detect a cold front in the 300° 
- 330° sector ~ 50" from the core, inside the region of 
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Table 2. Cold fronts candidates for the four clusters plotted in Fig. [3] We report the cluster name, the position angle (measured 
anticlockwise from East) of the candidate cold front, the different radial range considered, the slopes of the power laws that we 
find in the given radial range and the Aa values for the discontinuities in the surface brightness profiles. In the last column, a flag 
indicates which features are classified as discontinuities. 



cluster 


position angle 
(deg) 


radial range 
(arcsec) 


a 


ranges 


for Act 


Aa 


discont. 


Centaurus (internal) 


[120,150] 


10 


-50 


0.41 












Centaurus (external) 


[120,150] 


100 


-240 


0.84 












Centaurus (possib disc) 


[120,150] 


50- 


100 


2.59 


<*[50 


-100] 


- a[io-50] 


2.19 


V 


Centaurus (possib disc) 


[120,150] 


50- 


100 


2.59 


&[50- 


100] " 


' a [ 100-240] 


1.75 


V 


Centaurus (internal) 


[30,60] 


10- 


170 


0.88 












Centaurus (external) 


[30,60] 


210 


-600 


1.16 












Centaurus (possib disc) 


[30,60] 


170 


-210 


2.52 


»[170- 


210] " 


" ffpiO-600] 


1.64 


V 


Centaurus (possib disc) 


[30,60] 


170 


-210 


2.52 


«[170 


-210] 


~ C[10-170] 


1.36 


V 


A262 (internal) 


[-45,0] 


10 


-38 


0.84 












A262 (external) 


[-45,0] 


70- 


145 


0.72 












A262 (possib disc) 
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the weak shock wh ich is at ~ 200" dMcNamara et al.1120051 
iNulsen et aDl2005h . Similarly to A85, the sector of the cold 
front is narrow (30°). The jump is located at the bending 
of the south radio lobe (visib le in the Hydra A radio maps 
at 1.4 GHz; see for example lLane et ai1l2004h towards east, 
near the S W cavity. In that region the temperature map shows 
several cold blobs and one of this produces the temperature 
rise coincident with the surface brightness steepening but no 
clear front in the temperature map is present. The structure 
of Hydra A is very complex with a str ong interaction be- 
tween the radio lobes and the ICM gas (INulsen et al.ll2002t 
iMcNamara et alj2000l: see also lMcNamara & Nulsenl2 007). 
The cl uster exhibits a n umber of cavities in the central re- 
gions (I Wise et al.ll2.007h . The discontinuity we find is proba- 
bly a result of such a complicated morphology and likely it is 
not a cold front. Even if this region satisfies all the required 
conditions we consider this as an unclassified case. 

Since the existence of a cold front in A2052, A4059 and 
Hydra A cannot be definitively established, we exclude these 
systems from the sample. 

5. Cold fronts occurrence. 



5.1. Cold fronts occurrence: a general view 

The exclusion of three unclassified objects (namely A2052, 
A4059, Hydra A) reduces the sample to 42 objects of which 19 
host a cold front, corresponding to a fraction of 0.45. Note that 
the cold front in the A85 subcluster (dubbed A85*) is included. 
The list of the detected cold fronts, with their main properties, is 
filed in Table [3] 

Some clusters (e.g. Centaurus, A496, Perseus, A262) host 
more than one cold front. This phenomenon is not rare in cool 
core clusters. Chandra found multiple cold fro nts in several 
system s such as A 2 204, A 2029, Ophiuchus dSanders et all 
2005; Clarke et all l2004t IXscasibar & Markevitchl 12006: 
Markevitch & Vikhlinin 2007). The presence of multiple cold 
fronts in such clusters is likely related t o the origin and develop- 
ment of this phenomenon in cool cores t Ascasibar & Markevitchl 



l2006LlMarkevitch & Vikhlininll2007l) . We do not detect any cold 
front in some objects (i.e. A2204, A2029) , wh i ch are well- 
know n cold front systems (San ders et alj 120051 IClarke et alj 
120041) . For these clusters, the cold front is located in very central 
regions (14" and 30" from the center for A2029 and A2204, 
respectively). Therefore the discontinuity is well-resolved by 
Chandra, but it is hard to detect with XMM-Newton. 

A last comment concerns the bias in the measure of the 
occurrence of cold fronts due to instrumental and observa- 
tional limits: projection effects induce a smoothing on the sur- 
face brightness and temperature profiles and can hide a non- 
negligible fraction of cold fronts. Projection also completely 
hides cold fronts having an inclination larger than about 30° with 
respect to the plane of the sky. In addition, resolution limitations 
prevent the detection of cold fronts lying in the very central re- 
gions or cold fronts in distant clusters. Moreover, our detection 
algorithm may fail to detect some cold fronts with angular exten- 
sion < 30 °. All these effects significantly reduce the capability 
of detecting cold fronts. The frequency we find in our sample is 
therefore a lower limit of the true occurrence. 

5.2. Cold fronts occurrence: relation with redshift 

In our sample, no cold fronts are detected in systems at redshifts 
greater than about 0.075. We already remarked that A2204 and 
A2029 have been classified as cold front clusters from Chandra 
data analysis but we fail in detecting their discontinuities be- 
cause they lie at small distances ( < 30") from the X-ray peak, 
under the XMM-Newton resolution. For both A2204 and A2029 
Z > 0.075. This suggests that the lack of detection of cold fronts 
at the highest redshifts of the sample is most likely related to 
a resolution limit rather than to a real evolution. This effect is 
clearly shown in Fig. [4] where we plot the distances (in arcsec) 
from the cluster center of all the cold fronts detected in our sam- 
ple as a function of the redshift of the systems they belong to. 
Red points label merging clusters (see Table [3] and 35.3. Il l and 
black points label the remaining systems. Dashed dotted lines 
plot fixed physical distances (20, 50, 80, 150 kpc) at the various 
redshifts. From this figure, it is evident that cold fronts lying at 
~ 20-80 kpc from the cluster center are observed only in nearby 
systems (z ~ 0.05). Moving towards higher redshifts where these 
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Table 3. List of the cold fronts detected in the sample. For each cluster we report the center (RA and Dec) used to build surface 
brightness profiles, the azimuthal and radial positions of the cold front and the mean value of Aa. Bold faced fonts mark clusters 
having a merging cold front while italic fonts mark clusters whose merger geometry is not clear and the origin of the cold front is 
not as obvious (see 35.3b . 



Cluster name 


center 


position angle 


jump radius 


Ac 






(deg) 


(arc sec) 


(mean value) 


Centaurus 


12:48:49.173 -41:18:45.65 


[-60 , 60] 


170-210 


1.60 


Centaurus 


12:48:49.173 -41:18:45.65 


[90 , 210] 


50 - 100 


1.76 


A262 


01:52:46.117 36:09:05.79 


[-110, 0] 


38-70 


1.06 


A262 


01:52:46.117 36:09:05.79 


[30, 135] 


40-50 


1.40 


Perseus 


03:19:48 41:30:40 


[-60, 0] 


250 - 300 


1.72 


Perseus 


03:19:48 41:30:40 


[30, 150] 


140 - 200 


1.80 


A2199 


16:28:38.193 39:33:02.70 


[-75, -25] 


25-30 


0.50 


A496 


04:33:38.067 -13:15:40.91 


[-120, -30] 


40-55 


0.79 


A496 


04:33:38.067 -13:15:40.91 


[-120, -75] 


180 - 280 


1.33 


A496 


04:33:38.067 -13:15:40.91 


[30, 120] 


80 - 100 


1.86 


2A0335+096 


03:38:40.879 09:58:01.20 


[-120, -30] 


50-70 


1.16 


A 1644 


12:57:12.231 -17:24:32.67 


[-180, -90] 


20-35 


1.58 


A3558 


13:27:56.989 -31:29:50.00 


[-30, 120] 


80 - 120 


1.04 


A1795 


13:48:52.879 +26:35:27.80 


[-130 ,-80] 


60-70 


1.14 


A2065 


15:22:29.455 +27:42:23.81 


[-150, -90] 


80 - 100 


1.26 


A576 


07:21:30.495 +55:45:45.32 


[-120, -60] 


80 - 100 


0.72 


A3562 


13:33:36.766 -31:40:20.45 


[-150, -60] 


60 - 100 


1.06 


A1367 


11:44:53.5 19:44:19.12 


[-120 -60] 


350 ( 70 from the peak) 


0.84 


A754 


09:09:20.098 -09:40:52.22 


[120, 240] 


80 - 150 


1.14 


A85* 


00:41:42.733 -09:26:33.10 


[0, 90] 


25-80 


1.40 


A3667 


20:12:41.653 -56:50:52.94 


[-180, -90] 


250 - 280 


3.68 


A2319 


19:21:11.097 +43:56:08.00 


[-180, -60] 


100 150 ( 160 from the peak) 


1.85 


A2256 


17:02:33.009 +78:38:23.59 


[-135,-90] 


50 - 75 ( 100 from the peak) 


2.00 


A3266 


04:31:13.951 -61:27:26.41 


[-90, -30] 


60 100 


0.87 



physical distances progressively fall below the resolution limit 
(30", red solid line in the figure) cold fronts cannot be detected 
anymore. For z > 0.05, we have found only cold fronts at dis- 
tances > 80-100 kpc from the cluster center, with a prominent 
presence of merging systems whose cold fronts are generally lo- 
cated at large distances from the core (see 35.3.11 ). 

On the basis of the analysis of Fig. [4] we decided to apply 
a further selection on our sample, namely we cut the maximum 
redshift at z = 0.075 (i.e. the redshift where we stop detecting 
cold fronts). The resulting sample is reduced to 32 objects with 
a cold front occurrence of 59%. We note that the sample may 
be biased against clusters having cold fronts at small distances 
from the center, inducing an underestimation of the cold front 
occurrence we measure. However, in the case of cold fronts lying 
at distances r i 40 kpc the sample can be considered, to a first 
approximation unaffected by a redshift bias. 

5.3. Occurrence and origin of cold fronts. 

In this section, we investigate what discriminates clusters with- 
out cold fronts from clusters hosting one or more. 



5.3.1. Merger cold fronts 

We start by focussing our attention on merger cold fronts. Some 
of our systems are well known merging clusters. The mor- 
phology of these systems is generally complex and no unique 
center can be identified in their surf ace brightness maps. In 
many of these systems (e.g. A3667 Vikhlinin et al.l I2001 alibi 



IVikhUnin & M arkevitch 20ol lBriel et al.ll2 004) the merger pro- 
cess is occurring close to the plane of the sky and the geom- 
etry of the event is clear. The origin and the evolution of cold 
fronts in these systems is most likely related to the merger pro- 



cess. More specifically, the motion of a cold dense core of a 
subsystem which moves in the atmosphere of the main cluster 
during a merger event induces the formation of a cold front fea- 
ture. Typically, the subcluster is stripped of its outermost gas 
and the ram pressure exerted on the surviving dense cloud by 
the less dense surrounding gas produces the contact disconti- 
nuity between the two s ubsystems, generating the cold front 
dMarkevitch & Vikhlininll2007l) . 

In other objects, such as A3562 and A576, where the X-ray 
merger geometry is not as clear, the nature of the cold fronts 
we observe is not as obvious. A3562 is a cluster lying in the 
core of the Shapley supercluster, one of the largest mass con- 
centrations in t he local Universe . The p resence of a radio halo 
in this cluster (Giacintucci et al. 2005) provides an indication 
of a merger activity, since radio halos have been found only in 
interacting systems . Evidences of inter action come also from 
Bepp o-SAX data (iBardeiii et al.l [2002h . Using XMM-Newton 
data, [Finoguenov et al.l d2004l) suggest that the SC 1329-313 
group southwest of A3562 has passed to the north of A3562 and 
the cluster core is likely oscillating in response to the passage 
of the group. A576 is another peculiar syste m whose discontinu- 
ity ha s been observed also with Chandra by Kempn er & Davidl 
(120041) who propose that the core of the cluster is the remnant 
of a mergi n g sub cluster. This picture was also suggested by 
Mohr etall d 19961) from an analysis of the galaxy population. 
Recently, Dup ke etaLl ([2007) found that the system is consistent 
with a line-of-sight merger. According to this picture, the cold 
front we find in A576 is likely a merger cold front. 

The clusters of our sample hosting merger cold fronts are 
listed in the last part of Table [3] and marked with a bold-faced 
font. We include in this class A85*, the A85 subclump falling on 
the main structure. A576 and A3562, which are merging clusters 
where the cold front origin is not as readily associated to the 
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Fig. 4. Distances from the cluster center of the cold fronts de- 
tected in the sample plotted as a function of the redshift of their 
hosting systems. Red points label merging clusters (see Tabled 
and black points label the remaining systems. We omitted in this 
figure A85* since the cold front lies in a subclump. Dot dashed 
lines plot fixed physical distances at the various redshifts. Red 
solid line marks the XMM-Newton resolution limit at 30". 



Fig. 5. Scaled entropy profiles for the subsample described in 
95.41 Each curve is a locally weighted fit to the data points (see 
text for details) to reduce the scatter. Red solid lines are the pro- 
files of clusters hosting cold fronts, while green solid lines are 
the profiles of clusters without cold fronts. The red dot dashed 
curve is the profile for A3558. Cygnus A is omitted (see text for 
details). 



merger event, are placed in a separate category and marked with 
an italic font in Table [3] 

5.4. Non-merger cold fronts: the entropy profile 

In this subsection we investigate what discriminates clusters 
without cold fronts from clusters hosting at least one, once we 
exclude the clusters having a merger cold front (see above § 
15.3. U and clusters where the X-ray merger geometry is not clear. 
We focus on the remaining subsample (23 clusters) which in- 
cludes both clusters undergoing a merger event which is not ly- 
ing in the plane of the sky and clusters which do not present any 
sign of merging processes. In this subsample only 10 out of the 
23 clusters host a cold front (main properties are listed in the first 
part of Table[3]). We try to understand what determines the pres- 
ence of cold fronts in these systems, studying the radial entropy 
profile for each of these clusters. As is conventional in X-ray 
astronomy, we quantify the entropy using the adiabatic constant 
K = kTn~ 2 ^ (T and n e are the gas tem perature and dens ity re- 
spectively and A: is a constant) following lVoit et all ([2005). 

The specific entropy s is related to K through the relation 
s oc \nK. We will refer to K as "entropy" throughout the paper. 
To obtain the entropy profiles, we derive the radial profiles of the 
electron density n e and the temperature T, by deprojecting the 
observed surface brightness and temperature, under the assump- 
tion of spherical symmetry. The projected temperature and the 
surface brightness have been derived throug h a spectral analysis 
of the clusters using concentric annuli (see IRossetti & Molendi 
1201 Ol for details). To perform depro j ection , we adopted the pro- 
cedure described in iGhizzardiet ail (120041) . 



In Fig. [5] we plot all the derived entropy profiles. Radii are 
scaled to r^o, the radius within which the mean density is 180 
times the critical densit)Q. The values of r^o have been de- 
rived using its relationsh ip with th e cluster mean temperatur e 
as in lArnaud et al.ld2005l) (see also iLeccardi & Molendll l2008). 
The entropy is scaled using the empirical ent ropy scaling law 
K oc h(z) 4/3 T° ln 65 h 2 (z) = Q m (l + z) 3 + Q A dPratt et alJl200a 
iPonman et alj2003h : T\o is the mean tempe rature of each cluster 
in units of 10 keV, as in iPratt et alj ((2006). The curves plotted 
in Fig. [5] are a ctually a locally weighted re gression (LOWESS 
regression, see Sa nderson et alJuOOo, 12005) in the log-log space 
to reduce the scatter and provide a better view of the profiles be- 
havior. Clusters hosting a cold front (hereafter CF clusters) are 
denoted with red lines while clusters without cold fronts (here- 
after NCF) are denoted with green lines. We marked as particular 
case A3558 (red dot dashed curve). The behavior of this cluster 
will be discussed in Sec. [6] 

The cluster Cygnus A has been discarded here because of 
the strong contamination by the central AGN: the presence of 
the two hot spots invalidates the assumption of spherical symme- 
try and does not allow to deproject the temperature and surface 
brightness radial profiles. 

Fig. [5] shows that all the profiles have a similar trend at 
large radii (r > 0.08riso)- Moving towards the innermost re- 
gions, the profiles spread out and we observe a large scatter. 
More precisely, the NCF clusters (green solid curves) typically 
have central entropies higher than CF clusters (red solid curves). 



1 M m = mp c (z)(4n/3)r] m , where p c (z) = A 2 (z)3# 2 /8;rG and 
h\z) = n m (l + zf + fi A 
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Fig. 6. Thin lines are the mean scaled entropy profiles for clus- 
ters hosting cold fronts (red area) and clusters without cold fronts 
(green area). The shaded areas represent the standard deviation 
from the mean profiles. The thick blue line is the power law ob- 
tained fitting all the profiles in the radial range [0.08-0.3]ri8o- 
Blue points are the distances of the detected cold fronts (exclud- 
ing A3558) from the cluster center. 



Moreover, systems hosting cold fronts seem to have a steeper 
profile than clusters without cold fronts. 

In Fig. [6] we plot the two averaged profiles (thin solid curves) 
for CF and NCF clusters. The shaded areas represent the stan- 
dard deviation from the mean profiles. As in Fig. [5] color codes 
label CF clusters (red area) and NCF clusters (green area). We 
find that the mean profiles are similar at large radii. Fitting all 
the entropy profiles with a power law in the radial range [0.08 
- 0. 31nsn, we f ind a slope a = 0.95 ± 0.01. For compari- 
son^ 



Pratt 



et al] d20 

iPratt & Arnaudl 00051) find a = 0.94 + 0.14 and lPiffaretti et alJ 
(120051) find a = 0.95 + 0.02. Restricting to [0. 1 - 0.3]nso we find 
a slightly steeper power law with a = 1.08 ± .02 in accordance 
with th e theoret ical value of 1.1 predicted by |Tp zzi & No rmanl 
(120011) (see also lVoit & Ponmanll2003l:lBorgani et al.ll2002l) . 

Moving towards the innermost regions, for r H 0.08ri8o, the 
two mean profiles decouple. The NCF mean cluster profile ex- 
hibits a central entropy excess with respect to the outer power 
law model. The slope in the [0.01 - 0.08]rigo range is a — 
0.64+0.01 significantly lower than the outer power law slope. On 
the contrary, the CF cluster profiles become steeper in the same 
radial range, reaching lower central values. Fitting with a power 
law the CF mean entropy profile in the [0.01 - 0.08]ri8o range, 
we find a slope of 1 .22 + 0.0 1 , significantly higher than the exter- 
nal power law slope. The mean entropy values in the innermost 
bin (r = 0.005ri 80 ) are 94.5 ± 5.5 keV cm 2 and 10.8 + 1.6 keV 
cm 2 for NCF and CF clusters respectively. Excluding the partic- 
ular case A3558 (see $6]l does not significantly change results on 
best fit values. 



2006) found a slope a = 1.14 + 0.06, while 



Since the steepness of the entrop y profile is an indicato r 
of the presence of a "cool core" (e.g. Cavagn olo et alJ l2009). 
clusters with a steep entropy profile likely feature also a tem- 
perature decrement and a brightness excess in their internal re- 
gions. Therefore, one may argue that we do not detect cold fronts 
in objects with a fiat entropy profile just because their surface 
brightness is lower than that of clusters where we do detect cold 
fronts. However, as it can be seen from the profiles reported in 
Appendix, we detect cold fronts in the centers of cool core ob- 
jects (e. g. A262) where the SB is high, but also in the outer 
regions of merging clusters where the SB is much lower (e.g. 
A3667, A85*). In Fig. [3] we show an example of a cluster where 
we do not detect CF, AWM7, which has a SB comparable to 
those of the clusters where we do detect cold fronts. This is the 
case for almost all the clusters of our remaining sample where 
we do not detect cold fronts. 



6. Discussion 

The origin of cold fronts in clusters manifestly undergoing a 
merger event can be related to the motion of a dense cold cloud 
of gas within the atmosphere of another subcluster. Conversely, 
the presence (or the absence) of these features in the subsam- 
ple of non-merging clusters and clusters where the merger is 
not close to the plane of the sky is not clearly understood. Fig. 
[6] provides some hints to help understand what determines the 
presence of cold fronts in these systems. The general picture 
emerging here is that the entropy profile discriminates among 
the two classes (CF and NCF) of clusters. While at large radii 
the (scaled) entropy profiles of these clusters are very similar, in 
the innermost regions (r £ 0.08ri8o) their behaviors differ. 

Our finding of a steep entropy gradient in CF cluster is in 
agreement with theoretic al expectations. Indeed, simulations by 
Ascasibar & Markevitchl d2006l) show that cold fronts can rise 
and develop in the cores of clusters if the entropy sharply de- 
creases towards the center (as typically occurs in the center of 
cool core clusters). According to these simulations, cold fronts 
develop as a consequence of minor merger events; during its pas- 
sage near the center of a cluster, a merging subclump induces 
some disturbance on the low entropy gas of the core and dis- 
places it from the center of the potential well. If the entropy 
profile is steep, the cool gas starts sinking towards the mini- 
mum of the gravitational potential, a sloshing mechanism sets 
in and cold fronts arise. If the entropy profile is not steep, the 
entropy contrast is insufficient for the cool gas to flow back and 
for the sloshing mechanism to set in. In agreement with this pic- 
ture, we find that cold fronts form only in regions where the en- 
tropy profile sharply decreases. In Fig. |6j we plot (blue filled 
circles) the cold fronts positions measured for the non-merging 
clusters of our sample. Excluding the outermost cold front of 
A496 (this cluster hosts three cold fronts) which lies at a dis- 
tance of ~ 0.08ri8o from the peak, all the cold fronts we detect 
lie at distances smaller than ~ 0.05ri8o, where entropy profiles 
steepen, and greater than ~ O.Olriso where, in m any systems, 
Chandra detects a flattening (Donahue et al. 2003). 

The large majority of clusters of the final subsample obey the 
general rule that cold fronts are hosted by systems with a steep 
entropy profile in their centers. However, as already pointed out, 
A3558 is a peculiar case: although its entropy profile is similar to 
the NCF clusters profiles, it hosts a cold front. Some comments 
are needed to understand why this outlier does not follow the 
general behavior. The cold front for A3558 (blue square in Fig. 
13 is located at larger distance (r ~ 0.05ri8o) with respect to 
all the other cold fronts we detect, in a region where a weak 
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Fig. 7. Mean scaled entropy profiles as in Fig. [6] Red dot-dashed 
curve is the profile of A3558. The blue filled square is the A3558 
cold front distance from the cluster center. 



entropy gradient, not as sharp as in the other CF clusters profiles, 
is present. 

This cluster lies at the center of the Shapley supercluster and 
its special behavior of this cluster is likely related to the unique 
environment in which it is embedded. To understand the rea- 
son wh y cold fronts can arise in such a system, we refer once 
more to Ascasibar & M arkevitchl (120061) simulations. When the 
entropy profile of the cluster does not sharply decrease in the 
center, the central cold gas is easily pushed away from the dark 
matter peak, at the merging subclump passage. Accordingly, the 
cold front emerges at a large distance from the core. However, 
there is no entropy contrast to trigger the sloshing mechanism 
and this cold front will not develop further. Cold fronts rising in 
these syste ms are short-lived and therefor e rare phenomena (see 
Fig. 12 in lAscasibar & Markevitchll2006l) . A3558 is embedded 
in a very unrelaxed environment where merging events are fre- 
quent, and therefore the probability to form (and to observe) such 
fronts is higher. Alternatively, this cold front might be a merger 
cold front that we failed to recognize due to the fact that A3558 
cannot be classified easily as a merging cluster. As discussed in 
iRossetti et alJ d2007l) . it presents some features similar to those 
of cool core clusters and other properties that are more common 
in merging clusters. 

One of the main findings of our paper is that we detect 
at least one cold fro nt in all steep entropy gradient clusters in 
the final subsample. Ascasibar & Markevitch (2006) show that 
once the sloshing mechanism sets in cold fronts can be recog- 
nized in all the projection planes, even if they are more promi- 
nent on the merger plane (see Fig. 19 of their paper). However, 
the limited resolution of our instruments allow us to recognize 
only the most apparent brightness discontinuity. Indeed, we have 
performed some simulations of cold fronts projection with the 
XMM-Newton PSF and we found that cold fronts can only be 
observed if they lie within some 30° of the plane of the sky. This 



means that our 100% detection rate implies that most steep en- 
tropy clusters must host more than one "prominent" cold front. 
This abundance of cold fronts suggests that, whatever the trig- 
gering mechanism might be, it must have a high occurrence rate. 
Since the prominent cold fronts that we can detect are located 
on the merger plane, the detection of one or more cold fronts in 
all our steep entropy systems seems to indicate that a sizeable 
fraction of them are currently experiencing more than one minor 
merger. Assuming that, crudely speaking, cold fronts are visible 
for a timescale of about 3 Gyr (this is the case for the dark mat- 
ter + gas simulation in lAscasibar & Markev itch 2006, while for 
the dark matter only this timescale is longer), our cold front de- 
tection rate translates into a minimum merger frequency of 1/3 
merger event per halo per Gyr. If we further assume a minimum 
mass ratio of 1/10 we can compare our rate with rates expected 
from cosmologi cal simulations. Using Fig. 8 in a recent paper by 
Fakh oun "&Mal J20081) we find a merger rate of ~ 0.2 merger 
per halo per Gyr for mass ratio larger that 1/10. This is some- 
what smaller than the minimum rate implied by our observed 
cold front rates however, given the numerous simplifications we 
have applied in our calculation, we deem it to be in acceptable 
agreement. 

Gas sloshing may provide an important contrib ution to the 
coolin g-heating problem in cool core clusters (IZuHone et al.l 
120091) . The sloshing gas typically moves at sub-(or trans-)sonic 
velocities carrying a kinetic energy comparable to the thermal 
energy but the dissipation of this kinet ic energy to thermal en - 
ergy is too slow compared to cooling (Markevit ch et al.ll2001l) . 
However the sloshing mechanism also brings the outer high 
entropy gas into the core, mixing it with the cooling gas and 
resulting in a heat inflow which can prevent the formation of 
a "co oling flow" for periods of time 1-3 Gyr dZuHone et al.l 
2009). If subcluster encounters are frequent enough, as it is 
suggested by our high detection rate, the sloshing mechanism 
can efficiently offset cooling. Intriguingly the sloshing mech- 
anism operates preferentially in steep entropy profile clusters, 
i.e. precisely those which require heating to offset the cooling. 
With the coming into operation of the first space-borne micro- 
calorimeter, quite likel y the one onboard the ASTRO-H mission 
(Taka hashi et al.ll2008l) . it will be possible to investigate gas mo- 
tions in the direction of the line of sight, i.e. orthogonally with 
respect to that of the plane of the sky sampled with cold fronts. 
The combination of the two informations will afford a reliable 
estimate of the motions of the ICM in clusters core and estimate 
their role in offsetting cooling. 



7. Summary 

We have performed a systematic search of cold fronts using 
XMM-Newton data for a sa mple of 45 objec ts extracted from 
the B55 flux limited sample dEdee et alJll990h . 
The main results of our work are the following: 

- Excluding three unclassified cases, we find that 19 clusters 
out of 42 host at least one cold front. 

- We do not detect any cold front in systems having redshift 
greater than about 0.075. This is most likely related to XMM- 
Newton resolution limit. By cutting our sample at z = 0.075, 
we restrict our sample to 32 objects with a cold front occur- 
rence of 59% . 

- Cold fronts are easily detected in systems that are manifestly 
undergoing a merger event in (or close to) the plane of the 
sky. 
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- Out of the 23 clusters of the remaining subsample (systems 
undergoing a merger event which is not lying in the plane of 
the sky and non-merging clusters) 10 objects exhibit a cold 
front. For this final subsample, the entropy profile of systems 
hosting cold fronts is found to be steeper than that of clusters 
without them. The difference is observed at radii smaller than 
about 0.08ri8o where all our cold fronts are found. 

- Our findings are in agreement with simulation based predic- 
tions. As shown by lAsca sibar & Markevitch] d2006l) an en- 
tropy gradient is a necessary ingredient to trigger gas slosh- 
ing. 

- Since projection effects highly limit the capability of detect- 
ing cold fronts, the finding that all the clusters with a steep 
entropy profile host a cold front implies that most clusters 
with a steep entropy profile must have more than one cold 
front. 

- Under the assumption that cold fronts in cool core clus- 
ters are triggered by minor mergers, we estimate a mini- 
mum of 1/3 events per halo per Gyr, which is somewhat 
larger than that exp ected from cosmological simulations 
(Fak houri&Mall2008l) . 

- Gas sloshing may provide an important contribution to the 
cooling-heating problem in cool core clusters. A robust as- 
sessment of the gas motions associated to the sloshing phe- 
nomenon! will become possible with the coming into opera- 
tion of the first space borne microcalorimeter. 
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Appendix A: Surface brightness and temperature 
profiles 

In this Appendix we report the EPIC flux images 
of the clusters of the sample where we detected 
cold fronts (Full resolution maps are omitted in the 
astro-ph version of the paper and are available at: 
http://www.iasf-milano.inaf.it/~simona/pub/coldfronts/ghizzardi.pdf). 
The figures show the flux images in the 0.4-2 keV band, the 
black arcs indicate the position of the cold fronts and the "X" 
symbol the selected center for the extraction of the profiles. In 
Figures IA.UA.6I we show the surface brightness and tempera- 
ture profiles across the discontinuities in a representative sectors 
for all the cold fronts reported in Table [3] For each plot the 
discontinuity region is marked with vertical red dashed lines. 
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Fig. A.l. Surface brightness and temperature profiles for the cold fronts in Centaurus and A262. 
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Fig. A.2. Surface brightness and temperature profiles for the cold fronts in Perseus, A2199 and A496. 
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Fig. A.3. Surface brightness and temperature profiles for the cold fronts in A496, 2A0335 and A1644. 
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Fig. A.4. Surface brightness and temperature profiles for the cold fronts in A3558, A1795, A2065 and A576. 



16 



Simona Ghizzardi et al.: Cold fronts in galaxy clusters 



4 I 

. I 

• I 



A3562 (-150' ,-120°) 1 



10 100 
r (arsec) 



A3562 (-150° ,-130") 



200 300 
r (arsec) 





r (arsec) 



200 250 300 350 400 450 500 
r (arsec) 



A754 (120° ,150°) 



r (arsec) 



I i 



50 



A75^- (120° ,150°) 



100 150 200 250 300 
r (arsec) 



A85 (0' 



10 

r (arsec) 



A85 (0° ,90°) 



40 

r (arsec) 



Fig. A.5. Surface brightness and temperature profiles for the cold fronts in A3562, A1367, A754 and A85*. 
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Fig. A.6. Surface brightness and temperature profiles for the cold fronts in A3667, A2319, A2256 and A3266. 



